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Figure 1. Examples of the types of polygons studied in this paper. 



1. Introduction 

A well-known long standing problem in combinatorics and statistical mechanics is to 
find the generating function for self-avoiding polygons (or walks) on a two-dimensional 
lattice, enumerated by perimeter. Recently, we have gained a greater understanding 
of the difficulty of this problem, as Rechnitzer jl] has proved that the (anisotropic) 
generating function for square lattice self-avoiding polygons is not different iably finite 
[2], as had been conjectured earlier on numerical grounds [3]. That is to say, it 
cannot be expressed as the solution of an ordinary differential equation with polynomial 
coefficients. There are many simplifications of this problem that are solvable jl], but 
all the simpler models impose an effective directedness or equivalent constraint that 
reduces the problem, in essence, to a one-dimensional problem. 

A staircase polygon can be viewed as the intersection of two directed walks starting 
at the origin, moving only to the right or up and terminating once the walks join at a 
vertex. It is well-known that the generating function for staircase polygons is 



\ I - 2x - y/1 - Ax ^ , 2 _ a 

P(x) = 2 K ( ^ ' 

where the connective constant /x = 4 and the critical exponent a = 3/2. Punctured 
staircase polygons jH] are staircase polygons with internal holes which are also staircase 
polygons (the polygons are mutually- as well as self-avoiding). In [5 j it was proved that 
the connective constant /j of /c-punctured polygons (polygons with k holes) is the same as 
the connective constant of unpunctured polygons. Numerical evidence clearly indicated 
that the critical exponent a increased by 3/2 per puncture. The closely related model 
of punctured discs was considered in 0|. Punctured discs are counted by area and in 
this case it was proved that the critical exponent increases by 1 per puncture. Here 
we will study only the case with a single hole (see figure QJ, and we will refer to these 
objects as punctured staircase polygons. The perimeter length of staircase polygons is 
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even and thus the total perimeter (the outer perimeter plus the perimeter of the hole) 
is also even. We denote by p n the number of punctured staircase polygons of perimeter 
2n. The results of [5J imply that the half-perimeter generating function has a simple 
pole at x = x c = = 1/4, though the analysis in [3] clearly indicated that the critical 
behaviour is more complicated than a simple algebraic singularity 

Recently we found that the perimeter generating function of three-choice polygons 
can be expressed as the solution of an 8th order linear ODE [7] . Similarly, in this paper 
we report on work which has led to an exact Fuchsian linear differential equation of 
order 8 apparently satisfied by the perimeter generating function, V{x) = J2 n >oP nXn ^ 
for punctured staircase polygons (that is, V(x) is one of the solutions of the ODE, 
expanded around the origin). The first few terms in the generating function are 

V(x) = x 8 + 12x 9 + 94x 10 + 604X 11 + 3463x 12 + • • • . 

Our analysis of the ODE shows that the dominant singular behaviour is 

V(x) ~ A W + B(x) + C(x)bgtl-4*) + D(x){1 + 4x)13/2> (1) 

So in the notation used above, the generating function has a dominant singularity at 
x — x c — l/jj, with exponent a = 3. This result confirms exactly the conjecture for 
the critical exponent [5] in the case of a single puncture and the quite complicated 
corrections at the critical point explains why the analysis in j^J, based on a relatively 
short series, was so difficult. 

It is also possible to express the generating function V(x) as a sum of 4 x 4 Gessel- 
Viennot determinants [Hj. This is clear from figure El where the enumeration of the 
lattice paths between the dotted lines is just the classical problem of 4 vicious walkers, 
and these must be joined to two vicious walkers to the left, and to two vicious walkers 
to the right. Then one must sum over different possible geometries. However the fact 
that the generating function is so expressible implies that it is differentiably finite fSj. 

Unfortunately we cannot readily bound the size of the underlying ODE, otherwise 
we could use this observation to provide a proof of our results. As it is, we 
originally generated the counts of punctured staircase polygons up to perimeter 502 
(251 coefficients), and found what we believe to be the underlying ODE experimentally 
from the first 195 coefficients. The ODE then correctly predicts the next 56 coefficients. 
While the possibility that the underlying ODE is not the correct one is extraordinarily 
small, our procedure still does not constitute a proof of course. We have since extended 
the count beyond perimeter 800 and still all coefficients are predicted by our ODE. 



2. Computer enumeration 

The algorithm we use to count the number punctured staircase polygons is a modified 
version of the algorithm of Conway et al [10J for the enumeration of imperfect staircase 
polygons. The two problems are very similar and consequently there are only minor 
differences between the algorithms. A detailed description of the algorithm we used 
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Figure 2. The decomposition of a punctured staircase polygon into a sequence of 
2-4-2 vicious walkers, each expressible as a Gessel-Viennot determinant 

to count imperfect staircase polygons can be found in jjj. The algorithm is based on 
transfer matrix techniques. This entails bisecting the polygons by a line (as illustrated in 
figure and enumerating the number of polygons by moving the line 'forward' one step 
at a time. Punctured staircase polygons start out as ordinary staircase polygons and 
the line bisects the polygon at two edges. Then at some vertex two additional directed 
walks (sharing the same starting point) are inserted between the two original walks. The 
line will thus intersect these polygon configurations at four edges (see figure HJ). The 
only difference between the algorithm in |7] and the one used for this paper is in how 
the four directed walks intersected by the line are connected in order to produce a valid 
polygon. To produce a punctured staircase polygon we first connect the two innermost 
walks and then the two outermost walks are connected. Imperfect staircase polygons on 
the other hand are produced by connecting the first walk with the second walk and the 
third walks with the fourth walk. The updating rules used to count imperfect staircase 
polygons are given in [7j and are easily amended to count punctured staircase polygons 
bearing in mind the different 'closing' criteria outlined above. 

We calculated the number of punctured staircase polygons up to perimeter 502. The 
integer coefficients become very large so the calculation was performed using modular 
arithmetic ^1] . This involves performing the calculation modulo various prime numbers 
Pi and then reconstructing the full integer coefficients at the end. We used primes of the 
form pi = 2 30 — rj, where are small positive integers, less than 1000, chosen so that pi 
is prime, and Pi ^ pj unless i = j. 17 primes were needed to represent the coefficients 
correctly. The calculation for each prime used about 200Mb of memory and about 8 
minutes of CPU time on a 2.8 GHz Xeon processor. Naturally we could have carried 
the calculation much further (and we have since done this) but as we shall demonstrate 
in the next section this number of coefficients more than sufficed to identify an exact 
differential equation satisfied by V(x). 
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3. The Fuchsian differential equations 

In recent papers Zenine et al [T2"| IT3| ITij obtained the linear differential equations 
whose solutions give the 3- and 4-particle contributions x an d to the Ising model 
susceptibility. In [Zj we used their method to find a linear differential equation for three- 
choice polygons and in this paper we extend this work further to find a linear differential 
equation which has as a solution the generating function V(x) for punctured staircase 
polygons. We briefly outline the method here. Starting from a (long) series expansion 
for the function V(x) we look for a linear differential equation of order m of the form 



such that V(x) is a solution to this homogeneous linear differential equation, where the 
Pfc(x) are polynomials. In order to make it as simple as possible we start by searching for 
a Fuchsian equation. Such equations have only regular singular points. There are 
several reasons for searching for a Fuchsian equation, rather than a more general D-finite 
equation. Computationally the Fuchsian assumption simplifies the search for a solution. 
^From the general theory of Fuchsian equations it follows that the degree of Pk{x) 
is at most n — m + k where n is the degree of P m (x). To simplify matters further (reduce 
the order of the unknown polynomials) it is advantageous to explicitly assume that 
the origin and x = x c = 1/4 are regular singular points and set Pk(x) = Qk(x)S(x) k , 
where S(x) = x(l — 4x). Thus when searching for a solution of Fuchsian type there are 
only two parameters, namely the order m of the differential equation and the degree 
q m of the polynomial Q m (x). One may also argue, less precisely, that for "sensible" 
combinatorial models one would expect Fuchsian equations, as irregular singular points 
are characterized by explosive, super-exponential behaviour. Such behaviour is not 
normally characteristic of combinatorial problems arising from statistical mechanics. 
The point at infinity may be an exception to this somewhat imprecise observation. 

We then search systematically for solutions by varying m and q m . In this way we 
first found a solution with m = 10 and q m = 11, which required the determination of 
L = 195 unknown coefficients. We have 251 terms in the half-perimeter series and thus 
have 56 additional terms with which to check the correctness of our solution. Having 
found this solution we then turned the ODE into a recurrence relation and used this to 
generate more series terms in order to search for a lower order Fuchsian equation. The 
lowest order equation we found was eighth order and with q m = 27, which requires the 
determination of L = 294 unknown coefficients. Thus from our original 251 term series 
we could not have found this 8 th order solution since we did not have enough terms 
to determine all the unknown coefficients in the ODE. This raises the question as to 
whether perhaps there is an ODE of lower order than 8 that generates the coefficients? 
The short answer to this is no. Further study of our differential operator revealed that it 
can be factorised. In fact we found a factorization into three first-order linear operators, 
a second order and a third order. The generating function is a solution of the 8 th order 




(2) 



k=0 
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Table 1. Critical exponents for the regular singular points of the Fuchsian differential 
equation satisfied by V{x). 



Singularity Exponents 



x = 


-1, 0, 0, 0, 1, 2, 3, 8 


x = 1/4 


-1, -1/2, -1/2, 1/2, 1, 3/2, 2, 3 


x = -1/4 


0, 1, 2, 3, 4, 5, 6, 13/2 


x = ±i/2 


0, 1, 2, 3, 4, 5, 6, 13/2 


1 + x + 7x 2 = 


0, 1, 2, 2, 3, 4, 5, 6 


1/x = 


-2, -3/2, -1, -1, -1/2, 1/2, 3/2, 5/2 


Q 8 {x) = 


0, 1, 2, 3, 4, 5, 6, 8 



operator, not of any of the smaller factors. 

The (half)-perimeter generating function V(x) for punctured staircase polygons is 
a solution to the linear differential equation of order 8 

8 d fc 

k=0 

with 



P 8 (x) = x 4 (l - 4x) 8 (l + 4x)(l + 4x 2 )(l + x + 7x 2 )Q 8 {x). 

P 7 {x) = x 3 (l - 4:x) 7 Q 7 (x), P 6 {x) = 2x 2 (l - 4:xfQ 6 (x), 

P 5 (x) = 6x(l - AxfQ^x), P A {x) = 120(1 - 4x) 4 Q 4 (V), (4) 

P 3 (x) = 120(1 - Ax) 3 Q 3 {x), P 2 {x) = 720(1 - Ax) 2 Q 2 {x), 

P x {x) = 720(1 - 4a;)Q 1 (a;), P (x) = 2880Q {x), 

where Qs(x), Q 7 (x), . . ., Qo(x), are polynomials of degree 22, 28, 29, 30, 31, 31, 31, 31, 
and 31, respectively. The polynomials are listed in |Appendix A| 

The singular points of the differential equation are given by the roots of Ps{x). One 
can easily check that all the singularities (including x = oo) are regular singular points 
so equation (jSJ) is indeed of the Fuchsian type. It is thus possible using the method 
of Frobenius to obtain from the indicial equation the critical exponents at the singular 
points. These are listed in Tabled 

We shall now consider the local solutions to the differential equation around each 
singularity. Recall that in general it is known EH] that if the indicial equation yields 
k critical exponents which differ by an integer, then the local solutions may contain 
logarithmic terms up to log fe_1 . However, for the Fuchsian equation @ only multiple 
roots of the indicial equation give rise to logarithmic terms in the local solution around 
a given singularity, so that a root of multiplicity k gives rise to logarithmic terms up to 
log fc_1 . In particular this means that near any of the 22 roots of Qs(x) the local solutions 
have no logarithmic terms and the solutions are thus analytic since all the exponents 
are positive integers. The roots of Q$ are thus apparent singularities PHI of the 
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Fuchsian equation (J3J). We briefly mention that as in our earlier study [7] we can find a 
solution of order 14 of the same form as (jUJ) but with Qu(x) being just a constant. So at 
this order none of the roots of Qs{x) appear. Clearly any real singularity of the system 
cannot be made to vanish and so we conclude that the 22 roots of Q 8 must indeed be 
apparent singularities. 

Assuming that only repeated roots give rise to log-terms, and thus that a sequence 
of positive integers give rise to analytic terms, then near the physical critical point 
x = x c = 1/4 = l//iwe expect the singular behaviour 

A(x) B(x) + C(x) log(l-4s) 
/W ~(l-4x) + y/T=Ax~ ' [) 

where A(x), B(x) and C(x) are analytic in a neighbourhood of x c . Note that the terms 
associated with the exponents 1/2 and 3/2 become part of the analytic correction to 
the (1 — 4x) _1//2 term. Near the singularity on the negative x-axis, x = x_ = -1/4 we 
expect the singular behaviour 

V(x) ~ D(x)(l + 4x) 13/2 , (6) 

where again D(x) is analytic near x_. We expect similar behaviour near the pair of 
singularities x = ±i/2, and finally at the roots of 1 + x + 7x 2 we expect the behaviour 
E(x)(l + x + 7x 2 ) 2 log(l + x + 7x 2 ). 

We can simplify the 8 th order differential operator found above. We first found 
three very simple solutions of the ODE, each corresponding to an order one differential 
operator, 

= 1 - Ax, 

and 

„ . N 9 - 34x + 14x 2 
F 3 (x) = . . 

The existense of these three linearly independent solutions implies that we can find three 
first order operators, which we denote by l\ , with i = 1,2,3, such that the original 8'th 
order differential operator can be written as = L^L[ L 2 L% , where is a fifth 
order differential operator. We further found that this latter operator is decomposable 
as = L^'L^ 2 '. This then allows one, in principle, to write down the form of the 
8x8 matrix representing the differential Galois group of L^, in an appropriate global 
solution basis. To determine the asymptotics one would need to calculate non-local 
connection matrices between solutions at different points. This is a huge task for such 
a large differential operator. Instead, we have developed a numerical technique that 
avoids all these difficulties, which we describe below. 

To standardise our asymptotic analysis, we assume that the critical point is at 1. 
The growth constant of punctured staircase staircase polygons is 4, so we normalise the 
series by considering the new series with coefficients r n , defined by r n = p n+8 /A n . Thus 
the generating function we study is TZ(y) = J2 n >o r nV n = 1 + 3y + 5.875y 2 + ■ ■ -. Using 



Punctured staircase polygons 



8 



the recurrence relations for p n (derived from the ODE) it is easy and fast to generate 
many more terms r n . We generated the first 100000 terms and saved them as floats 
with 500 digit accuracy (this calculation took less than 15 minutes). ^From equations 
(0) and © it follows that the asymptotic form of the coefficients is 

— = E (|+^^+(-ir [jm)) + o(a-»). ( 7, 

Any contributions from the other singularities are exponentially suppressed since their 
norm (in the scaled variable y = x/4) exceeds 1. 

Estimates for the amplitudes were obtained by fitting r n to the form given above 
using an increasing number of amplitudes. 'Experimentally' we find we need about 
the same total number of terms at x c and — x c . So in the fits we used the terms with 
amplitudes di, bi, and q, i = 0, . . . , K and di, i = 0, . . . , 3K. Going only to K with 
the di amplitudes results in much poorer convergence and going beyond 3K leads to no 
improvement. For a given K we thus have to estimate 6K + 4 unknown amplitudes. So 
we use the last 6K + 4 terms r n with n ranging from 100000 to 100000 — 6K — 3 and 
solve the resulting 6K + 4 system of linear equations. We can also add extra terms to 
the asymptotic form and check what happens to the amplitudes of the new terms. If 
these amplitudes are very small it is highly likely that the terms are not truly present 
(if the calculation could be done exactly these amplitudes would be zero). 

Doing this we found that all the amplitudes di appear to be zero except that 
d Q = 1024, e.g., with K = 20 we find that the estimates for the amplitudes d\ < 10 -70 , 
<i2 < 10 -60 , ds < 10~ 50 , etc. So in all likelihood the amplitudes di = for i > 0. This 
then leads us to the refined asymptotic form 

M = r „ = 1 024 (l + -i=W^^ + (- 

V V i>0 v 

In fits to this form we then used the terms with amplitudes bi, and Cj, i = 0, . . . , K 
and di, i = 0, . . . , 2K. For a given K we thus have to estimate 4K + 3 unknown 
amplitudes. We find that the amplitude estimates are fairly accurate up to around 
2K/3. We observed this by doing the calculation with K = 30 and K = 40 and then 
looking at the difference in the amplitude estimates. For b and Cq the difference is less 
than 10 -120 , while for do the difference is less than 10~ 116 . Each time we increase the 
amplitude index by 1 we lose around six significant digits in accuracy. With i = 18 the 
differences are respectively around 10~ 14 and 10~ n . 

^From our very long series it is possible to obtain accurate numerical estimates 
of many of the amplitudes bi, Ci, and di, with a precision of more than 100 digits 
for the dominant amplitudes, shrinking to around 10 digits for the the case when 
i = 18 (actually we could probably have pushed this further but there would be 
little point). In this way we found that b = -|^|, h = 4v g 3/2 , b 2 = 192 8 S /2 , 
c = 1.55210340048879105374 . . . and d = J, A = -fp, d 2 = ™, d 3 = 
di = ^2Q4g^3/^ ; though we have not been able to identify c . These amplitudes are 
known to at least 100 digits accuracy. 
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The excellent convergence is solid evidence (though naturally not a proof) that 
the assumptions leading to equation (J7|) are correct. Further evidence was obtained as 
follows: We have already argued that the terms of the form 1/ri 1 , i > are absent. We 
found similar results if we added terms like \ogn/n\ log 2 n/n^ 2 or additional logn terms 
at y — — 1. So this fitting procedure provides convincing evidence that the asymptotic 
form (jHJ), and thus the assumptions leading to this formula, are correct. 

4. Conclusion and Outlook 

We have developed an improved algorithm for enumerating punctured staircase 
polygons. The extended series, coupled with a search program that assumes the solution 
is a Fuchsian ODE, enabled us to discover the underlying ODE, which was of 10 th order. 
We did this without using 56 of the coefficients that we had generated. That is to say, 
56 of the known coefficients were unused, and so their values provided a check on the 
solution. This leads us to believe that we have found the correct ODE, as it reproduces 
the known, unused coefficients. Subsequently we have extended this check to more than 
200 unused coefficients. Further refinement allowed us to find an 8 th order ODE. 

A numerical technique we have developed specifically for such problems then allowed 
us to find accurate numerical estimates for the amplitudes of the first several terms in 
the asymptotic form of the coefficients, and then to conjecture their exact value. 

We have also initiated an investigation of the area generating function. We expect 
this to involve g-series, and thus far our investigations only lead us to believe that the 
area generating function A(q) is of the form 

A(q) = {G{q)+H{q)y/T^^)/[J (l,l,q) 3 \, 

where J (x,y,q) is a g-generalisation of the Bessel function, and occurs, for example, 
in the solution of the problem of staircase polygons enumerated by perimeter ^j. Here 
q = r) is the first zero of Jo(l, 1, q), and G and H are regular in the neighbourhood of 
q — rj. The coefficients thus behave asymptotically as 

a n = [q n ]A(q) ~ const.rf n n. 

In a subsequent publication we propose to investigate the area generating function more 
fully, and hopefully obtain more insight into the properties of the ODE we have found 
for the perimeter generating function. 

Furthermore in work with C. Richard ^7] we have conjectured the scaling function 
for punctured polygons with an arbitrary number of punctures. We briefly review the 
properties of the two- variable area-perimeter generating function for staircase polygons. 
Of special interest is the point (x c , 1) where two lines of singularities meet. The 
behaviour of the singular part of the generating function about (x c , 1) is expected 
to take the special form 

P(x, q) ~ P^\x, </) + (!- qfF((x c - x){l - q)~% (x, q) /, 



Punctured staircase polygons 



10 



where F(s) is a scaling function of combined argument s = (x c — x)(l — q)~^, commonly 
assumed to be regular at the origin, and 6 = 1/3 and = 2/3 are critical exponents. 
For staircase polygons, we have 

F( a ) = ~togAi((V2)*a). 

In ^7j we assumed that the limit distribution by area of staircase polygons with r 
punctures (of arbitrary size) is that of staircase polygons with r holes of unit area. 
From this and knowledge of F(s) we then obtained exact predictions for r punctured 
staircase polygons by taking the r-th derivative w.r.t q of P(x,q). We then study the 
area-moment generating function, Pk(x) = J2 mn n k p m , n x m , where p m ,n is the number of 
polygons with perimeter m and area n. In particular we find that the leading amplitudes 
of the perimeter generating function of the k-th area-moment are given by 

( P) = (k + r)\x r J k+r 
k r\x2 k+r T( lk+r ) 

Here f k + r are amplitudes occuring in the asymptotic expansion of P(x, q) (these are 
known exactly for punctured staircase polygons) while ^k+r = 3(A; + r)/2 — 1/2 are the 
critical exponents of the fcth area-moment of r punctured polygons. These predictions 
have been confirmed for once punctured staircase polygons to a very high level of 
accuracy for moments up to k = 10. The numerical analysis of the area-moments 
relied crucially on our knowledge of the singularity structure of the perimeter generating 
function as detailed in this paper. 



E-mail or WWW retrieval of series 



The series for the generating functions studied in this paper can be obtained via e-mail 
by sending a request to I.Jensen@ms.unimelb.edu.au or via the world wide web on the 



URL http://www.ms.unimelb.edu.au/~iwan/ by following the instructions. 
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Appendix A. Polynomials Q n (x) for punctured staircase polygons 

Q 8 (ar) = - 90720 + 1255590a; - 9538200a; 2 + 20394828a; 3 - 79106610a; 4 

+ 1223958687a; 5 - 2943232056a; 6 + 17470357067a; 7 - 189472079743a; 8 

+ 579172715823x 9 - 1746461498616x 10 + 8400325324610X 11 

- 1591154327260a; 12 - 111431714394808a; 13 + 315517552430480a; 14 



Punctured staircase polygons 

- 106489387477312a; 15 - 938487878760384a; 16 + 1628517397980288a; 

- 2394531569420032a; 18 + 2966185168205312a; 19 

- 170238270849024a; 20 - 699187344629760a; 21 + 295462090506240a; 2S 



Q 7 (x) = - 1360800 + 23565780a; - 167569290a; 2 + 478254996a; 3 + 641052858a; 4 

+ 14810951034a; 5 - 47034372339a; 6 - 166933659974a; 7 - 2552936187594a; : 

+ 64479 11404224a; 9 + 14253364474478a; 10 + 86598771199392a; 11 

+ 362131239586500a; 12 - 3860712252484892a; 13 + 8993313236994576a; 14 

- 31235880957264960a; 15 + 46429326957124912a; 16 

+ 155905775680790304a; 17 - 807736441103822976a; 18 
+ 1835072857042276096a; 19 - 1278888252797142528a; 20 

- 293981468599460352a; 21 + 14541716059525437440a; 22 

- 26481815895022608384a; 23 + 22483566008412450816a; 24 

- 35911819535956066304a; 25 + 3639680241277796352a; 26 
+ 7495959535363031040a; 27 - 3507725938490081280a; 28 



Qe(x) = - 1723680 + 69281730a; - 787195710a; 2 + 4886678970a; 3 - 10726639974a; 4 
+ 11830409583a; 5 - 401281487235a; 6 + 343905413598a; 7 
+ 1858137414650a; 8 + 44092692217413a; 9 - 36740412036168a; 10 

- 135298590380414a; 11 - 1279093006602396a; 12 - 10004750418032976a; 13 
+ 61536871579988144a; 14 - 216281351081049504a; 15 

+ 1050287576547538488a; 16 - 1795967175346626976a; 17 

- 2572736181692580960a; 18 + 18017037664470796032a; 19 

- 45232775265352713472a; 20 + 48709527110201501184a; 21 
+ 4770083118869915136a; 22 - 322327838255331590144a; 23 

+ 541571044899035842560a; 24 - 511926023257614434304a; 25 
+ 716375351150156644352a; 26 - 69659801950830723072a; 27 

- 136551990333116252160a; 28 + 60094625512245166080a; 29 



Q 5 (x) = 1965600 - 6539400a; - 358033410a; 2 + 4831433820a; 3 - 30915098190a; 4 
+ 60211846008a; 5 - 201764518161a; 6 + 2531858233470a; 7 
+ 1380416576424a; 8 - 20212314275250a; 9 - 61506470769366a; 10 
- 477804842150324a; 11 + 608746761166938a; 12 + 483723642457152a; 13 
+ 60127368616743592a; 14 - 185780400624937008a; 15 
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+ 1165835175099337288a; 16 - 7175943616536571776a; 17 
+ 13745698284061066112a; 18 + 4948349174336379840a; 19 

- 89453290124304769024a; 20 + 270104157697832561664a; 21 

- 356324521463829808128a; 22 - 41862184650482117632a; 23 

+ 1845216328946812827648a; 24 - 2906213125616330383360a; 25 
+ 2943265956913569742848a; 26 - 3723507915329643413504a; 27 
+ 405249143061461336064a; 28 + 618215144006850969600a; 29 

- 261821958729561538560a; 30 



Q 4 (x) = 241920 - 8017380a; + 88351704a; 2 - 590355612a; 3 + 2409400818a; 4 

- 8457027588a; 5 + 71232186468a; 6 - 288557341128a; 7 

+ 524905454055a; 8 - 5046532132734a; 9 + 28114089314043a; 10 

- 164508486596467a; 11 + 869331744354740a; 12 - 2401501341116904a; 
+ 12275987679372578a; 14 - 50846889626226508a; 15 

+ 46258831828476364a; 16 - 147764159295056304a; 17 
+ 1375769527659995736a; 18 - 2625251094439093408a; 19 

- 765792895039661984a; 20 + 22951686058011476032a; 21 

- 85054223999548283904a; 22 + 126294091912315062016a; 23 
+ 19381267403906712064a; 24 - 566287434634380073984a; 25 
+ 849895463062111623168a; 26 - 892557255237919469568a; 27 
+ 1043719341871898804224a; 28 - 142670999896790335488a; 29 

- 140350544778022354944a; 30 + 59234239904690995200a; 31 



Q 3 (x) = - 4596480 + 112443660a; - 1327020156a; 2 + 11580963786a; 3 - 76436209584a; 4 
+ 426159579924a; 5 - 2350462539072a; 6 + 11395385983233a; 7 

- 44136036344190a; 8 + 145288111685523a; 9 - 559910802106640a; 10 
+ 3013037795053530a; 11 - 13499762948930634a; 12 

+ 50096716464628528a; 13 - 217987216302493908a; 14 
+ 853439326193439492a; 15 - 2363497210984795232a; 16 
+ 5096223845046539304a; 17 - 8508469151526998016a; 18 
+ 9581930085552894304a; 19 - 10241374665198721536a; 20 

- 12641088914996048640a; 21 + 118651673978481267200a; 22 

- 208768950 136609496064a; 23 - 15400291418459486208a; 24 
+ 814317146169694152704a; 25 - 1202858442211165741056a; 26 
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+ 1271933402411862171648a; 27 - 1406355411740766470144ar 
+ 251165051564655771648a; 29 + 137326949251639934976a; 30 
- 61285928661166325760a; 31 



Q 2 {x) = 1209600 - 10784340a; + 25225200a; 2 - 192390408a; 3 + 1497608946a; 4 

- 3085618896a; 5 + 55270573062a; 6 - 674664767886a; 7 + 1891951243653a;' 
+ 6937954472784a; 9 - 19443421819978a; 10 - 252270853719194a; 11 

+ 1421753108033868a; 12 - 2280488850916676a; 13 - 1040351739238056a; 1/ 

- 1519080794794788a; 15 + 54144924827952720a; 16 

- 143110935850986376a; 17 - 63031554528921744a; 18 

+ 1125126938486807936a; 19 - 2675665192031509504a; 20 
+ 3361130538055156224a; 21 - 2669659667713374208a; 22 
+ 1996890960732463104a; 23 - 4866848788151009280a; 24 
+ 3555378162093901824a; 25 + 3193922372633202688a; 26 

- 2642707373157531648a; 27 + 2132642211038560256a; 28 

- 3311881541411143680a; 29 + 1596569887904366592a; 30 

- 264734033093591040a; 31 



Q^x) = - 725760 + 19969740a; - 254689092a; 2 + 2329185726a; 3 - 17948325636a; 4 
+ 118028863386a; 5 - 679983561900a; 6 + 3637871524611a; 7 

- 17150360490738a; 8 + 62088405193554a; 9 - 183555964459890a; 10 

+ 747009873725220a; 11 - 4106684548673028a; 12 + 18540613780587884a; 13 

- 67936944600058776a; 14 + 247341581626824360a; 15 

- 939866071520217104a; 16 + 3216462341735279616a; 17 

- 8789133587934808704a; 18 + 17976423995943224576a; 19 

- 26625353996773725696a; 20 + 29354499014436664320a; 21 

- 26197184327864145920a; 22 + 20118012206750361600a; 23 

- 11595016904008224768a; 24 - 12803308242930466816a; 25 
+ 49275320633035751424a; 26 - 49679788190366564352a; 27 
+ 31169615491025600512a; 28 - 29010025645678264320a; 29 
+ 12772559103234932736a; 30 - 2117872264748728320a; 31 



Qo(x) = Qi{x) 
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